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Abstract 

A consensus that questions the perfunctory use of the quantum adiabatic theorem has emerged 
since Marzhn and Sanders [Phys. Rev. Lett. 93, 160408 (2004)] showed the existence of an 
inconsistency in the apphcabihty of the theorem. Further analysis proved that the inconsistency 
Q_| may arise from the existence of resonant terms in the Hamiltonian, but recent work indicates that 

the debate about the full extent of the problem remains open. Here, we first show that key premises 
required in the standard demonstration of the theorem do not hold for a dual Hamiltonian involved 
in the Mar zlin- Sanders inconsistency. Also, we show that two simple conditions can identify systems 
for which the adiabatic approximation fails, in spite of satisfying traditional quantitative conditions 
that were believed to guarantee its validity. Finally, we prove that the inconsistency only arises 
for Hamiltonians that contain resonant terms whose amplitudes go asymptotically to zero. 

PACS numbers: 03.65. Ca, 03.65.Ta, 03.65.Vf 
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I. INTRODUCTION 



A folk quantum adiabatic theorem (QAT) establishes that a system initially described 
by an eigenstate associated with an isolated eigenvalue, of an instantaneous time-dependent 
Hamiltonian, H{0), will be at further times in the corresponding eigenstate of H{t), if H 
changes slowly enough. The rate at which the adiabatic limit is reached depends on the gap 
separating the eigenvalue of interest from the rest of the energy spectrum. 

The theorem was initially proven by Born and Fock |l| for bounded Hamiltonians with 
discrete energy levels and extended later by Kato who removed the assumption of spectral 
discreteness, provided the initial eigenstate corresponds to a discrete eigenvalue. On the 
other hand, Avron and Elgart [sl showed that an adiabatic theorem can be formulated for 
systems with no spectral gaps. In practice, when the Hamiltonian variation is not infinitely 
slow, the adiabatic theorem does not hold although, under certain conditions, the dynamics 
can be approximately adiabatic. Statements about these conditions constitute the so-called 



adiabatic approximation j4|. The present paper deals, somehow simultaneously, with the 
theorem and the approximation, but we emphasize that both concepts are not synonymous 
and they must be clearly distinguished. 

The publication in 2004 of a paper sl pointing out an inconsistency in the theorem has 
resulted in a new consensus among an increasing number of physicists, which indicates that 
"extreme caution should be used when interpreting results based on the standard application 



of the theorem" |6|]. The so-called Marzlin-Sanders inconsistency js] can be presented in the 
following simple way due to Tong et al. [3]: If a system Sa evolves adiabatically under the 
dynamics generated by a Hamiltonian Ha{t), a dual system S^, described by the Hamiltonian 
Hbit) = —Uait, Oy Hait)Uait, 0), whcre Ua(t, 0) is the exact evolution operator for system Sa, 
does not evolve adiabatically in general. However, both systems satisfy the same quantitative 
requirements 



« 1, m 7^ n , Vt , (1) 



Enit) — Emit) 

that traditionally were considered to guarantee that a given Hamiltonian (with eigenpairs 
En(t), \En(t))) evolves slowly enough. 

Tong et al. [7] attributed the inconsistency to the insufficiency of Eq. ([1]). Comparat 
js] argued that Eq. ([1]) is valid to detect adiabatic systems except when the Hamiltonian 



2 



contains oscillatory terms. MacKenzie et al. ^ \ emphasized the differences between adiabatic 
approximation and adiabatic theorem, and explained that the inadequacy of the quantitative 
conditions arises in situations where the approximation, but not the theorem, is invalid. Duki 

n 

et al. [9|] showed that the paradox results from the breakdown of the adiabatic approximation 



for the dual system. Similarly, Amin 10| showed that the inconsistency arises because the 



dual Hamiltonian contains resonant terms. Other authors [llj, 12| indicated that the dual 
Hamiltonian, if;,, present multiple time scales and for that reason adiabatic theorems do not 
apply to it. Finally, some works have denied the inconsistency [isl 14 1, 

The current consensus indicates that QAT has been rigourously proven beyond any doubt, 
and it is not in question. In fact, no problems have been reported so far in the derivations of 
the theorem. The dominant notion seems to be that the Marzlin-Sanders paradox is related 
to a failure or insufficiency of the conditions that assure that the adiabatic approximation 
is applicable to a particular system. FoUowingthis line of thought, large efforts have been 



dedicated to propose new adiabatic conditions [4 
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. Nonetheless, recent work 



21| indicates that the debate about the full extent of the inconsistency remains open. 



Here, by carefully examining the standard proof of QAT due to Messiah 22|, we show 
in Sec. |TT]that the theorem does not hold for the dual system because some key premises 
required in the demonstrations of the theorem are not satisfied by Hh. This confirms that 
only the adiabatic approximation is affected by the Marzlin-Sanders inconsistency, while 
the theorem is correct ^ . As a result of our analysis we show in Sec. IIIII that two simple 
conditions can be used to determine if a system for which Eq. ([T]) holds satisfies the adiabatic 
approximation. We illustrate our findings in Sec. HVlwith an example. In Sec. |V]we discuss 
the full extent of the inconsistency, showing that it only arises when the Hamiltonian contains 
resonant oscillatory terms whose amplitude goes to zero in the asymptotic limit. Finally, 
we give a summary and a general discussion in Sec. |VT1 



II. THE QUANTUM ADIABATIC THEOREM 



Premises 



The adiabatic theorem refers to the limit in which a time-dependent Hamiltonian Ha{t) 
varies infinitely slowly. The theorem gives approximate solutions to the time-dependent 
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Schrodinger equation in an interval < t < r, when r — )■ oo and the total change of Ha{t) 
is finite [2]. Using a scaled time variable s defined hy t = rs, < s < 1, the Schrodinger 
equation can be written 

= THa{s)<P{s) . (2) 

In the following, we assume that the energy spectrum is discrete at all times as the 
Marzlin-Sanders paradox arises in such a context. Hence, the following theorem 



Ua{s)P:XO) - P:{s)Uais) = 0(l/r), r ^ oo, Vn , (3) 

can be proven under the following additional premises 22|]: (pi) the eigenvalues E'^{s) are 
continuous functions of s, (p2) there are no eigenvalue crossings (i.e., E^{s) ^ E^{s), n m, 
s G [0, 1]), and (p3) the derivatives of the eigenprojectors, dP^{s)/ds and d'^P^{s)/ds'^ , are 
well defined and are piecewise continuous in the interval < s < 1. For simplicity Ha{s) is 
usually assumed to be independent of r but the arguments leading to the adiabatic theorem 
apply similarly to the case in which H„ depends explicitely on r, as far as premises (pl)-(p3) 
are uniformly fulfilled for r — )■ oo {2], |2^ . 



B. Sketch of the standard proof of the theorem 



Messiah's demonstration 



22| is based on the use of a virtual evolution operator, Ua{s) 



defined by Kato [2], that takes the system from a given eigenstate at the initial time to the 



corresponding eigenstate at future times. A purely geometric evolution is realized by 



24| 



UAis) = j2ms)){E:{o)\, (4) 

n 

where the phases of the instantaneous eigenfunctions, |£'"(s)), are such that the parallel- 
transport condition {E^{s)\E^{s)) = is satisfied. Operator Ua fulfills the so-called inter- 
twining property 

Ua{s)P:{0) = P:{s)Ua{s) , (5) 
and therefore, QAT, Eq. ([3]), can be proven by showing that 
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lim U\{s)Ua{s) = ^a{s) [1 + 0(l/r)] , (6) 

T— ^OO 

where the operator ^a{s) can be expanded in terms of the eigenprojectors at s = multiphed 
by phase factors 



J]exp \{-iT/n) rE:ia)da 
„ I Jo 



pm . (7) 



Messiah showed that the solution of the evolution equation for operator VIa = U\Ua goes 
to ^A in the limit r — )■ oo. For that, he made two unitary transformation over the evolution 
equation in the Schrodinger representation. This procedure leads to the following Volterra 
integral evolution equation for operator Wa = ^^^^a 

Wa{s) = 1 + (t/h) [ Ka{cT)WA{a)da , (8) 
Jo 

where the kernel, 

Kais) = ^\{s)U\{s)Ka{s)UA{s)^A{s) , (9) 

depends on the following operator 

Ka{s) = zhj2p:i^)p:i^), (10) 

n 

which is the generator of the virtual evolution given by Ua- The completion of the proof 
requires showing that the integral in the Volterra equation goes to zero when r — t- oo. If 
that is the case, Eq. ([8]) gives Wa 1, and from the definition of Qa results Qa — ^ ^a, and 
finally Ua Ua^a- 

An estimate of the integral in Eq. ([8]) can be obtained by substituting in Eq. ([9]) 
the standard expansion Ka = '^nmPni^)^aPm{0), and Eq. ([7]), which gives the following 
expression for the kernel 



Kais) 



ih exp 



ir/h) / {E^^{a)-E:^{a))da 



{El,{s)\E:{s))\E^M{Km- (11) 



Nondiagonal elements of the kernel contain exponential factors whose phases are never sta- 
tionary because they oscillate at frequencies that increase with r, as far as the differences 



E^{a) — E!^{a) are nonzero (i.e., if the eigenvalues are separated by gaps). Messiah showed 
that the integral in Eq. ([8]) can be expressed as a sum of two terms, which contain as a 
factor the operator 



FJs] 



Ka{(y)da 



(12) 



Thus, if Fa{s) goes to zero, the integral in the Volterra equation will go to zero too. Due 
to the presence of oscillatory factors in the kernel, operator Fa{s) approaches zero if the 
elements {E'^{s)\E'^{s)) are continuous functions of s, and their derivatives remain finite for 
all s |22|]. 

III. IDENTIFICATION OF SYSTEMS AFFECTED BY THE INCONSISTENCY 
A. The dual system 



Amin 



10| showed, by taking into account the adiabaticy of Sa, that the dual Hamiltonian 



Hb can be expanded as 



H,{s,T) = -Y.^n{s)\Em){Km 



ih \ ^ 

2^ exp 



[-ir/h) / {El{a) - El{a))da 
'o 



{E:^{s)\K{s))\E:,m{Ei 



(13) 



where the expression given in 



10| has been rewritten here as a function of scaled time. 



Similarly, the eigenprojectors, P;,, can be written as 



ih 



:-ZT/h) /Vn(^) - EUcr))da 
Jo 



x{Kis)ms))\Euo)){E:,m, (14) 

which shows that the Hamiltonian and the eigenprojectors for the dual system depend 
explicitely on r. This fact does not invalidate a priori the adiabatic theorem as far as 
premises (pl)-(p3) hold when r — t- oo. 
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Premises (pi) and (p2) depend only on the gaps between pairs of eigenvalues. Therefore, 
they are equally satisfied for Ha and Hi, since E^{s) = — However, premise (p3) does 
not hold for H^, since the operator- valued functions, dP^/ds and d'^P^/ds^, are not defined 
when r — t- oo, because the argument of the exponential in Eq. f[T^ goes to infinity. The 
explicit r-dependence of the Hamiltonian, when premise (p3) does not hold, prevents the 
adiabatic limit for Sb from being reached when r — )■ oo, unlike for Sa- 

The adiabatic theorem holds for the dual system if the integral in the Volterra equation, 
Eq. dH]), for Sb goes to zero in the r — ^ oo limit. As explained in Sec. IIIBI such an 
integral goes to zero whenever (i) the eigenvalues are separated by gaps and (ii) the elements 
(£'^(s)|-E'^(s)) are continuous functions of s that remain finite along with their derivatives 
for all s. Thus, the adiabatic theorem would be satisfied for the dual system if (ii) holds 
since (i) is automatically guaranteed because Hb{s) has the same spectral gaps as Ha{s), 
which is adiabatic by hypothesis. 

A wrong reasoning indicates that the integral in the Volterra equation should go to zero 
for system 5*;, as it does for Sa, since Eqs. (H1)- (|TT1) are apparently valid for Sb if subscript 
A{a) is changed to B{b). Then, the exact Qb would tend to an operator analogous to Eq. 
© and QAT would hold for Sb- From E^^{s) = -E^{s) and P^(0) = P^{0), this operator is 
= 

B. Properties satisfied by the dual system in the r — > oo limit 

The naivety of the previous argument is exposed by realizing that a virtual adiabatic 
operator for Sb is given, in general, as Ub = UIUa^- Note that a Volterra equation analogous 
to Eq. ([8]) is obtained only if Ub is a geometric evolution operator; that is, if it can be 
expanded like Eq. (jlj) in terms of eigenf unctions of Hb that satisfy {E'^{s)\E'^{s)) = 0. The 
following relationship between parallel-transport eigenfunctions of Sa and Sb, 

|E^(s,r))=exp 

allows us to expand the virtual evolution operator, Ub, in terms of eigenfunctions of Sa, 
which satisfy {E^{s)\E^{s)) = 0, as 



-tr/h) rE:^{a)da 
Jo 



ume:{s)) 



(15) 
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Ub{s,t) = J^exp 



-^T/h) / E:{a)da 



(16) 



On the other hand, the relation between the eigenprojectors of both systems 



p::{s,T) = ui{s)p:{s)u,{s), (17) 

gives Kb{s) = Ul{s)Ka{s)Ua{s), which can be substituted, along with $b(s) = $^(s), and 
Eq. f|T6l) . in the analogous equation to Eq. IQ appropriate for system Sb, to obtain 



Kb = ^lUlKbUB^B = ^A^\U\UaUlKaUaUlUA^A^ 



t 

A ' 



(18) 



where the s dependence has been omitted for simplicity. Finally, by using the equations, for 
Sb, analogous to Eqs. (jlj) and (ITOl) . Kb can be expanded as 



Kb{s) = zhJ2{E:n{s)ms))\E:,{0)){E^„ 
which does not contain oscillating functions. Hence, 



(19) 



Fb{s) = / Kb{a)do ^ 



r 



(20) 



which implies that the integral in the Volterra evolution equation does not go to zero. In 
other words. 



Kb{a)WB{(y)da ^ 



r 



(21) 



because 



{El{a)\E'Sa))da = r{E'^^{a)\K{a))da , 
Jo 

(22) 

due to the cancellation of the exponential factors. 

The adiabatic theorem does not hold for system Sb because the mechanism that makes the 
integral in Eq. (|8]) approach zero for Sa fails for Sb in spite of the fact that both systems have 
the same spectral gaps. The reason is that the oscillations of the elements {E'^{s)\E'^{s)) 
cancel the terms that arise from the gap condition. We emphasize that this does not imply 



lim 



exp 



{^r/h) / {El{a') - E^{a'))da' 
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an inconsistency in QAT since system Hh does not satisfy the hypotheses of the theorem. 
Here, we prove that systems for which the effects of the gap condition are canceled can be 
easily identified by simple conditions. Also, we will show that such conditions can be applied 
without modification to identify systems that do not satisfy the adiabatic approximation in 
spite of satisfying Eq. ([T]). 

By substituting the relationship between the eigenprojectors of systems Sa and Sb, Eq. 
f lT7|) . into the intertwining property, Eq. ([5]), which Ua satisfies due to the adiabaticity of 
Sa, we get 

[U1{s)Ua{s)] P„^(0) [UA{s)^Ua{s)] = P„^(., r) . (23) 
By taking the r — )■ oo limit, and substituting Eq. into the previous expression, we obtain 



limP„^(s,r) = P„^(0)[l + O(l/r)]^ Vn, s G [0, 1] , (24) 

which holds for systems affected by the Marzlin-Sanders inconsistency; that is, for systems 
that satisfiy Eq. ([1]) but do not evolve adiabatically. However, Eq. fl2^ . is ambiguous and 
therefore is not a sufficient condition. In effect, a system Sc with Hamiltonian —Hf, has the 
same eigenprojectors but, being the Heisenberg representation of Ha, is trivially adiabatic. 
This ambiguity can be eliminated by finding some other property that holds for Sb but not 
for Sc- The kernel for = —Hb is 

Kc{s) = ^\^\u\KaUA^A^A , (25) 

which shows that the phase cancellation that produces the inconsistency for system Sb does 
not occur for system Sc- The kernel expansion for Sc contains oscillatory factors that are 
related to those of Sa by 



exp 



Jo 



exp 



{Ki'y)\Ki'y)) , (26) 



(2^r/n) r{E:^{a') - E^^{a'))da' 
Jo 

where the eigenvectors l-E^), l-E'^), fulfill the parallel-transport condition. Since system Sa 
is adiabatic by hypothesis, the integral of the expression at the right-hand side of Eq. (!26|) 
goes to zero as 1/r, which implies, for r — )• oo. 



{E^M\K{'^))dcr = 0{l/ 



T 



(27) 



exp (zr/n) / {E^^ia') - E^S^'))da' 
Jo 

while, as shown in Eq. ( l22|) . the analogous integral for system Sb does not go to zero. 

In conclusion, if the eigenprojectors of a given Hamiltonian fulfill Eq. (IMll . the traditional 
adiabatic conditions, Eq. ([T]), are insufficient to determine the adiabaticity of the time 
evolution. In such a case, only if the integral in the left-hand side of Eq. f l22|) goes to zero, 
will the evolution be adiabatic. 



C. Properties satisfied by the dual system for finite r 

The adiabatic approximation, as mentioned in the Introduction, is the set of conditions 
under which the time evolution of a system that evolves during a finite time interval, is 
still approximately adiabatic. A system for which the adiabatic theorem holds when r — t- cxo 
could fail to evolve adiabatically for a certain set of finite values of the parameter r. However, 
if the adiabatic theorem does not hold when r oo, the evolution cannot be approximately 
adiabatic for any finite r. Since the dual system does not obey the adiabatic theorem it 
does not obey the adiabatic approximation either. However, the failure to obey premise 
(p3), which invalidates the theorem, refers to the strict adiabatic limit and cannot be used 
to determine whether the adiabatic approximation does not hold. Next, we show that 
there exist common conditions that apply to both the adiabatic theorem and the adiabatic 
approximation, from which we will obtain properties that can be used to determine the 
validity of Eq. 

Previous studies of the Marzlin-Sanders inconsistency have considere 
opment of new criteria of validity for the adiabatic approximation 0, [sl 
the following, we use directly the variable t for time instead of s, as it is traditional in the 
analysis of the adiabatic approximation. Appropriate criteria must identify unambiguously 
if a Hamiltonian that acts during a finite time interval varies slowly enough for the state 
vector at time t to be approximated by the eigenstate \En{t)) except for a phase factor, if 
the state vector at initial time was \En{Q)). 

In the r — )■ oo limit, conditions, Eqs. ( l22l) . and ( l24l) . are sufficient to indicate if a given 
Hamiltonian, iff,, does not satisfy the adiabatic theorem. These same conditions can identify 
systems that satisfy Eq. ([T]) but that, however, do not satisfy the adiabatic approximation. 



main 
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y the devel- 
, 15, 161. In 
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In effect, if Ha satisfies the adiabatic approximation, we get 

Ua{t) ^ UAmAit) . (28) 

Substituting this expression in Eq. (125]) . which approximately holds for finite but long 
enough r, we obtain 

Pnit) ^ P'M , Vt . (29) 

On the other hand, Eq. ( |22|) holds too, because the matrix elements {E^{t)\E^{t)) cancel 
the exponentials arising from the gap condition for finite r as in the r — t- oo case. Finally, 
if the adiabatic approximation is obeyed, the integral in Eq. (1271) will be small. 



IV. EXAMPLE 



The previous results can be illustrated by analyzing a model that has been extensively 
studied in relation to the Marzlin-Sanders inconsistency. This model consists of a spin-half 
particle in a rotating magnetic field whose Hamiltonian is j?! 

Hait) = —(X'o/2((j2; sin 6^ COS + CTy sin 6* sin + COS 6*) , (30) 

where uj and ujq are constants, while the cr-s are Pauli matrices. The instantaneous eigen- 
values are -E" (t) = LiJq/2 and E^it) = — a;o/2, while the instantaneous eigenvectors are |7| 



Ei{t))=\ M , (31) 




and 



-iujt/2 9 

mt))= I ' \ , (32) 

I e-*/2sinf ' 
which do not obey the parallel-transport phase condition. 

The Hamiltonian, Eq. (l30l) can be written as a function of scaled time s = 27rt/r = ut 

Ha{s) = —uo/2{ax sin 9 cos s + ay sin 9 sin s + az cos 9) . (33) 
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Eigenf unctions of Ha{s) that obey the parallel-transport condition are given, after trans- 
forming Eqs. ([3l]) and ([32]) to scaled time, by \En) = ^o^^^^\^^)<^''\En) ■ Thus, we get 



_g«s/2 ^,Qg . 



and 



-is/2 e 

Finally, note that the adiabatic limit for Ha{s) corresponds to r — > oo, which is equivalent 
to — )■ 0. Such a limit is not physically realizable, but that is not relevant for the illustration 
of the inconsistency. 



A. The w ^ limit 

We showed above that a dual Hamiltonian Hf){s,T) = —Ul{s,T)Ha{s)Ua{s,T) does not 
satisfy, in general, the adiabatic theorem. Let us study the case corresponding to Ha given 
by Eq. ( I33l) . Hamiltonian Hb can be written as 

Hbis,r) = -J2KPnis,r), (36) 

n 

where the instantaneous eigenprojectors are given by = UlP^Ua- Transforming the 
expression for Ua, given by Tong |7|, to scaled time, we get 



(cos ^ + i '^+^0_cose ^\ g-is/2 ^^^..^ gj^ ^g- 



I uiQ sin 6 gjj_|^ tJSg— js/2 

" ■■^-|iii^sin^e^^/2 (cos ^ - z^i^i^i^^ sin ^) 

u) 2u) \ 2u) uj 2uj / 



where w = a/wq + + 2wwo cos 6. On the other hand, the following matrix representation 
for can be constructed from the instantaneous eigenvector, Eqs. ([31 



sin^ I — sin I cos |e , 

A"(^)= I , % . I • (38) 

— sin I cos le*"* cos | 
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The matrix representation for the instantaneous eigenprojector, Pi{s), can be calculated 
from the previous two equations and has a complicated form. For example, one of its matrix 
elements is 



[P^{s,t)],2 = sin' ^—s 



Uq [U + Uq cos ( 



sin 6 sin' 



9 / Wo sin 6* ' 



6 9 
sm - cos - 
2 2 



/ u + uo COS 9\ .9 9 Wo sin 6* (cj + Wo cos 9) o 

+ sm - cos 7^ cos 

V w / 2 2 



w2 



ofuj\.9 9 
— cos — s sm - cos - 
V2w / 2 2 



. . f UJ \ f UJ 

+ « sm — s cos — s 
\2u J \2u 



uo . „ . 2^ , r,^ + '^ocos9 . 9 9 u;o 

— smt7sm — \- 2 sm - cos smt^cos - 

- 2 uJ 2 2 UJ 2 



(39) 



We showed in Eq. (fl^ that the eigenprojectors of the dual system contain terms that, 
although they oscillate infinitely fast, go to zero in the r — )■ oo limit. Apparently, this is not 
the case for [P^]i2- However, Eq. ([39]) can be written, by using well known trigonometric 
identities such as 



a;o\2 9 9 fu 
— sm - cos — h — 
u) 2 2 \uj 



■ 9 9 
sin - cos - 
2 2 



2 f W 



, , — \ ■ 9 9 „ . . 

— cos — s sm - cos — V2i sm 
' " ' 2 2 



V2w 

By taking into account that a; — implies cJ — )• wq, we get 









. 9 


9 


l2^^r°^i 






sm - 


cos - 




\2u, J 




2 


2 



(40) 



\nn[P',{s,T)]u 



9 9 
sm - cos - 
2 2 



. 2 / W \ f^\' ^- ■ I ^ \ UJ \ UJ 

1 — sm — s — + 2« sm — s cos — s — 
\2uj / \ujJ \2uj / \2uj / UJ 



= [Pmi2[l + 0{uj) + 0{uj')\ . (41) 
Following the same procedure for the other elements of the eigenprojector matrix we obtain 



limPi'(s,r) 



sm 



— sm 2 cos 2 



sm 2 cos 2 



cos 



l + 0(w)M(^.)+0(w^)N(ii. 



(42) 
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where M and N represent matrices whose elements are functions that oscillate infinitely fast. 
However, they are multiplied by terms that go to zero in the w — )■ limit. Equation f H2|) 
shows that -Pi (s) [and P2{s), which is not given here], in the w — )■ limit, tend to -Pi (0) [and 
P2^(0)]. Note also that the instantaneous eigenprojectors do not contain several independent 
time scales. The time dependence has the form ujs/{2uj), which suggests the use of a new 
scaled variable s' = ujs/{2u). Also, the eigenvalues are constant, for this particular example, 
and consequently the Hamiltonian can be written as a function of s' too. This shows that 
the adiabatic limit for Hf,, for which QAT would hold, requires u/ {2ijj) — )■ in addition to 
— > 0. This limit is physically meaningless. 

As explained above, Hamiltonian He = —Hb, has the same eigenprojectors as Hb but, is 
adiabatic. Both systems can be distinguished by evaluating 



r 




/ exp 


{i/uj) 1 


Jo 


Jo 



t/io) / {Ei{a') - E2{cT'))da' 



{E,{a)\E2ia))da 



(43) 



where En{s) = —E^{s) for Hb and -E„(s) = E'^{s) for He- From the relationship between 
the eigenfunctions of both systems and those of Ha we have 



{Ei{a)\E2i(T)) = exp 
and taking into account 



(z/u) {El{a') - E1{a'))da' 





{El{a)\E^,{a)) 



(44) 



Eq. ( I43l) becomes for Hb 



(45) 



exp 



[i/u) I {El{a') - E',{a'))da' 





{El{a)\El{a))da 



{El{a)\E^{a))da = ^(1 - e^^'=°^'') tan^ , 
which shows that Sb is not adiabatic unless ^ = 0, mod tt. However, for He-, we have 

{El{a)\E^,{a))da 



(46) 



exp 



[i/u) / {El{a') - E^,{a'))da' 
Jo 



exp 



{2t/io) I {E1{a') - E^{a'))da' 





{E'^iam{a))da 
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which goes to zero in the hmit w — )■ 0, indicating that the adiabatic theorem holds for system 
Sc, hke for system Sa, regardless of the values of the parameters Uq and 6. This example 
shows the validity of Eqs. ( |22|) . and (!24|) to identify if a Hamiltonian Hb is related to an 
adiabatic Hamiltonian Ha through Hb = —U\HaUa- 



B. The finite oo case 



System Sa satisfies the adiabatic approximation if uj is small compared to wo? which 
implies cJ ~ wq. If we employ this approximation in the expression for the matrix element 
of the eigenprojector Pi(t), which can be obtained writing Eq. (139|1 as a function of t, we 

get 



[Piit)] 



e e 

12 ~ — sm - cos - 
2 2 



2 f^o.\ . 2 f^o. 

sm^ I —tj + COS^ yt 



, . . , Wo \ /uJo.\ „ .06 
+1 sm — t cos — t X U = — sm - cos - 
2 / V 2 / 2 2 

(48) 



The other matrix elements can be simplified in the same way, so we finally obtain (t) 
Pi{0), which proves that Eq. fl29|) holds for this system. 

On the other hand, for finite u, i.e., for finite r, Eq. (H6|) becomes 



/ {E^{a)\E^{a))da = -(1 - e*'^*™^'') tan^ 
Jo 2 



(49) 



which indicates that system Sb satisfies the adiabatic approximation only for 6 = 0, mod vr, 
while the integral in Eq. (147|) is 



exp 



2i I (El {a') - E^{a'))da' 





{El{a)\E^,{a))d<T 



u sin 6[1 — exp (it (2(X'o + u) cos 6*))] 



(50) 



4(wo + 2wcos6') 

which is small for ujq » uj\ that is, system Sc satisfies the adiabatic approximation under 
the same conditions as system Sa, while 5*;, does not. 
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V. FULL EXTENT OF THE MARZLIN-SANDERS INCONSISTENCY 



Traditional adiabatic conditions, given by Eq. ([T]), hold for the dual Hamiltonian defined 
by Tong [3]. However, the time evolution driven by Hf, is not adiabatic. We have proven that 
the ultimate reason behind this inconsistency is the existence, in the Hamiltonian, of reso- 
nant terms that go asymptotically to zero. This makes the gap condition, usually invoked 
to guarantee adiabatic behavior in the r — )■ oo limit, become irrelevant. These peculiarities 
of Tong's dual Hamiltonian, Hh, are due to its special relationship with an adiabatic Hamil- 
tonian, Ha- In effect, both Hamiltonians are connected by a unitary transformation and a 
sign change. The transforming unitary operator is the exact time-evolution operator, Ua, for 
the adiabatic Hamiltonian. It seems pertinent to investigate the status of the inconsistency 
for Hamiltonians, H^, that are related to an adiabatic Hamiltonian, Ha, through a more 
general unitary transformation, with or without a sign change (i.e., = ±U^HaUx, where 

U,, ^ Ua). 



A. The family is generic 

At first sight it may seem that is affected by a fundamental restriction, since its 
eigenvalues must satisfy = ±i?^. However, we argue here that all possible cases of 
interest can be studied within this approach. Specifically, can contain generic oscillatory 
terms, as can be seen by turning upside down the previous argument relating and Ha- 
In other words, instead of arguing that generic oscillatory Hamiltonians can be obtained by 
unitary transformation of a particular adiabatic Hamiltonian, it can be argued that for a 
given Hamiltonian, H^, that by hypothesis is generic, appropriate unitary transformations 
lead to an adiabatic Hamiltonian. These unitary transformation are nonperturbative (i.e., 
they are not close to the identity). The resulting adiabatic Hamiltonian can be understood 



as an effective Hamiltonian that contains the effect of resonances 25 |. 

Therefore, given a Hamiltonian Hx{s,t) that is generic, in the sense that it may con- 
tain strong oscillatory terms (instead of the weak resonances that appear in the dual 
Hamiltonian used in the Marzlin-Sanders inconsistency), there exists an unitary oper- 
ator, Ux{s,t), that depends explicitely on r, such that the transformed Hamiltonian, 
±Ux{s,t)Hx{s,t)U^{s,t) = Ha{s), is adiabatic. Note that Ha is not necessarily the same 
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adiabatic Hamiltonian, Ha, studied in previous sections, since 7^ Ht,. Note also that the 
unitary operator Ux{s, r) is not the evolution operator for systems Sa- Thus, this procedure 
does not impose a priory any restriction on H^. However, the unitary operator, U^, is 
somehow restricted, because it must transform the generic into an adiabatic Hamilto- 
nian. This transformation can be very complicated but it is always possible for Hermitian 
matrices because they can be diagonalized. Thus, for a given H^, there is always a similarity 
transformation with an unitary operator, Ux = Ud, that converts into a diagonal matrix, 
D. The time evolution driven by H^, for an initial function can be expressed, in terms 



of Ud, as 



26| 



^h-^^'^{s) = TD{s)^is) - thUUs, t)^^^^^{s) , (51) 

where \l/ = Uj^^. When contains strong oscillatory terms the derivative of Ud will be 
large and the time evolution is not adiabatic. However, the evolution driven by Hamiltonian 
D is 

in^^(s) = tD{s)^{s) , (52) 

OS 

which is adiabatic since D does not contain off-diagonal couplings. This proves that a generic 
Hx can be related by an unitary transformation, driven by Ux, to an adiabatic Hamiltonian 
Ha- Strictly, the proof implies that, at least for Ux = Ud, and Ha = D, the transformation 
exists. However, it is clear that other adiabatic Hamiltonians can be generated from the 
same Hx by using less restrictive Ux operators. 



B. General relations between Hx and Ha 

A generic unitary operator Ux is always the evolution operator for an unknown Hamilto- 
nian H 7^ Ha 

th^^^^^ = TH{s,T)Ux{s,T). (53) 

Parallel-transport eigenfunctions, \E^), for Hx, can be obtained from eigenfunctions with 
arbitrary phase, \E^), as follows 



17 



\K) = exp 



- / (K\K)da 



\K)- 



Since we can chose \E^) = we get, by using Eq. (ES]), 



(54) 



\Kis,T))=exp 
Thus, for n ^ m, 



-ir/h) / {E^S^)\H{a,T)\El{a))da 
Jo 



Ul{s,T)\El{s)) . (55) 



{E:,{s,r)\El{s,r)) = (itT/h){E:;,is)\H{s,T)ms)) + {E:,{s)ms)) 



X exp 



[tr/h) 1^ (^{EM\Hia,T)\E:,ia))-{E:{a)\H{a,T)\EM))da 



(56) 



The relation between geometric evolution operators for Sr,. and Sa contains, in general, 
an additional phase factor, $, i.e., Ux = U}.Ua^- Thus, the kernel, Eq. ([9]), for is 

= $t,$^^i?7,ir,f/]^A<f $x , (57) 

where $x = if Ha = -UxH^Ul or $x = $a if Ha = Un^HJJl. Operator $x gets 
canceled, in the kernel, if $ = or $ = ^\ depending on the sign of the transformation 
connecting H^ and Ha- In both cases, the resulting Ux is not a geometric operator, except 
in special cases, because the parallel-transport phase fixing condition is not satisfied for 
eigenfunctions \E^) = f/j!$^|_E^) or \E^) = U}.^\\E'^). Hence, the kernel for a generic 
contains oscillatory terms. Also, Eq. (!57|) depends on the sign of the transformation that 
relates H^ to Ha'. 

K^{s,r) = -Tj2{Kis)\H{s,T) + {th/T){d/ds)\KmE:nm{Km 



xexp 



[tr/h) [{El{a)\H{a,r)±Ha{a)\El{a)) - {E^^{a)\H{a,T) ± Ha{a)\El{a))) da 

(58) 



We will use these relations to study the adiabaticity of Sx based on the characteristics of 
Ux- Then, we will compare the results with the left-hand side of Eq. ([1]) for S^- For that, it 
is necessary to write Eq. (|56|) as a function of usual time t instead of s, which gives 
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E^{t)-E^it) 



{i/h){E-^{t)\H{t)\E-{t)) + (E« 



(59) 



E'^it) - E-^(t) 

We present, in the following, heuristic arguments instead of rigurous theorems. Thus, 
this analysis should be seen only as a sketch for future developments of the present work, 
whose core is the material contained in Sec. IIIII 



C. Status of the inconsistency for various cases 

1. Hx does not contain oscillatory terms 

As Hx has, by hypothesis, spectral gaps, if it does not contain oscillatory terms, premises 
pl)-(p3) in Sec. |TT]hold. Then, in the r — t- oo limit, the evolution is adiabatic . Comparat 
8| showed that in this case Eq. ([1]) holds. Thus, no inconsistency a la Marlinz-Sanders 
takes place. In effect, adiabaticity for implies that 



Jo 



(60) 



which, in turn, implies that the kernel cannot contain terms that depend linearly on r. This 
means that matrix elements {E^{s)\H{s,t)\E^{s)) ~ 0. Thus, H must be close to Ha or it 
must be a constant with small time-dependent terms that vary slowly. The first case leads 
to a contradiction, since then the arguments of the exponentials get canceled and the system 
will not be adiabatic. Therefore, H must be given by a nearly constant term. Hence, Eq. 
( |59|) gives 



E^{t)-E^{t) 
and no inconsistency results. 



{Eum^it)) 



Em - Em 



« 1 ,yn,m;\/t , 



(61) 



2. Hx contains strong resonant oscillatory terms 

Previous studies j^, [lo| have related the Marzlin-Sanders inconsistency to the existence 
of resonances but have not discriminated between weak and strong perturbations. These 
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previous works did not consider either the r — )■ oo hmit. Thus, the full extent of the incon- 
sistency and its relation to the standard proof of the adiabatic theorem was not analyzed 
before. 

t i s well known that the existence of exact resonances precludes adiabatic behavior 
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27|. Thus, we have 



K^{(7)da ^0{1/t) . (62) 

On the other hand, Ux contains strong oscillatory nondiagonal terms, since it relates an 
adiabatic Hamiltonian to a Hamiltonian with strong oscillatory terms. Thus, H = ihUxUl 
will contain strong oscillatory terms too, and will be very different from Ha- In spite of the 
fact that the arguments of the exponentials in the kernel, Eq. (1581) . are not canceled, the 
kernel integral does not go to zero because the matrix elements {E^{s)\H{s,t)\E^{s)) are 
nonmonotonic functions that depend on r. 

The question at stake here is if Eq. ([1]) holds under these circumstances. If it does, there 
will be an inconsistency. Note the difference with Tong's dual system, for which = Ua, 
H = Ha, and the matrix elements {E'^{t)\H\E'^{t)) = 0. Instead, when the resonances are 
strong these matrix elements will not be small. Then the matrix elements {E!^{t)\E!^{t)) 
will not be small either. Therefore, the condition, Eq. ([1]), will not be satisfied, and no 
inconsistency results. 



3. Hx contains strong non-resonant oscillatory terms 

Nonresonant oscillatory terms can be safely averaged out, and therefore the exact time 
evolution is driven by the average Hamiltonian, which does not contain oscillating terms 
anymore. Thus, in the r — )■ oo limit the average Hamiltonian will be adiabatic. Numerous 
cases exist in the bibliography showing the adequacy of this approach. A well-known example 
is the time evolution of the molecular alignment that takes place when a molecule interacts 
via its polarizability with a strong nonresonant laser pulse. This interaction depends on 
the square of the electric field. The time evolution of the rotational wave functions can 
be faithfully studied by taking into account only the laser envelope after averaging out the 



rapid oscillatory terms cos^((X't), where u is the nonresonant laser frequency 28|]. Thus, if 



the laser pulse is long enough the evolution is adiabatic, and Sec. IV C II applies. 
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Summarizing, the Marzlin-Sanders inconsistency takes place for a system Sx only if the 
kernel of the Volterra evolution equation does not oscillate at frequencies that increase 
with T and does not contain terms that vary nonmonotonically. The kernel, Eq. 
fulfills these two conditions if if ± Ha = 0. The case H = Ha implies, from Eq. flS5]) . 
Ux = Ua, and the resulting Hx is precisely the dual Hamiltonian used by Tong [5] (i.e., 
Hx = Hb = —UlHaUa)- Contrarily, for the positive sign choice in the transformation that 
relates Hx and Ha-, the kernel satisfies the two conditions ii H = —Ha- In this case, Ux 
is the exact evolution operator f/^) for the reversed Hamiltonian —Ha- In rigour, the 
kernel fulfills the two cited conditions not only if if = iii^ but also if H admits a full 
asymptotic expansion, in powers of r~^, in which the zero-order term is given by ±iia 
[i.e., if H{s,r) = ±Ha{s) + J2ji'^^/'^y H'' i^)]^ so the properties derived in Sec. IIIII are 
approximately valid. 

VI. DISCUSSION 

The so-called Marzlin-Sanders inconsistency can be described as follows. If Eq. ([1]) is 
satisfied for a Hamiltonian Ha for which the adiabatic approximation holds, it will be satisfied 
too for the Hamiltonian Hf, = —UlHalla, where Ua is the exact evolution operator for Ha- 
However, in general, Hb does not satisfy the adiabatic approximation [7]. This implies that 
Eq. ([I]) is not a sufficient condition for a system to hold the adiabatic approximation. 

The full extent of the inconsistency has been much debated and doubts about the con- 
sistency of the own adiabatic theorem have been raised. Here, we have proven that these 
doubts cannot be sustained. Given a Hamiltonian Ha{s = t/r) for which an adiabatic theo- 
rem holds in the r — t- oo limit, the dual Hamiltonian Hb{s, r) = — f/](s, T)Ha{s)Ua{s, r) does 
not satisfy the theorem because Hb contains resonant terms that oscillate infinitely fast al- 
though their amplitude go asymptotically to zero. Due to these terms the dual Hamiltonian 
does not vary infinitely slowly when r — )■ oo. Also, the derivatives of the eigenprojectors 
are not defined. Thus, no inconsistency affects the quantum adiabatic theorem because the 
dual Hamiltonian does not satisfy key premises required in the proofs of the theorem. 

In Sec. Ill Al we showed that premises (pl)-(p3) were imposed to guarantee that some 
oscillatory integrals that arise in the evolution equations, due to the existence of spectral 
gaps, go to zero in the adiabatic limit. Therefore, they only make strict sense in such a limit 
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and cannot be used to determine the validity of the adiabatic approximation. However, the 
adiabatic approximation does not hold for the dual system because the same integrals that 
do not go to zero in the adiabatic limit are not small for finite r. On the other hand, the 
conditions Eq. ([T]) hold for the dual Hamiltonian because the matrix elements {E^{s)\E^{s)) , 
although they oscillate fast, have a small magnitude. 

Systems affected by the inconsistency can be easily identified because the instantaneous 
eigenprojectors of the dual Hamiltonian oscillate with negligible amplitude, with respect 
to a baseline, at frequencies resonant with the energy levels of the system. This baseline 
is defined by the eigenprojectors at the initial time. As a consequence, the instantaneous 
eigenprojectors change very little with time. However, the condition, Eq. (1241) . is neces- 
sary but not sufficient for the inconsistency to occur, because it is equally satisfied by the 
Hamiltonian —Hb, which is adiabatic. 

On the other hand, the presence of fast oscillations of very small amplitude does not 
invalidate the adiabatic approximation unless they cancel the oscillatory factors arising from 
the gap condition (i.e., unless they are resonant). Thus, resonant terms, even if very small, 
make the gap condition irrelevant. 

Also, we have proven, by using an approach different to that of Comparat j^, that 
the adiabatic approximation does not hold for Hamiltonians with strong oscillatory terms 
(resonant or nonresonant). In addition, our analysis indicates that, in this case, Eq. ([1]) 
is not satisfied. Hence, we conclude that Eq. ([T]) is a necessary and sufficient condition 
for a system to satisfy the adiabatic approximation unless the time-dependent part of the 
Hamiltonian contains resonant terms of very small magnitude. 

In cases for which an analytical expression is known for a particular Hamiltonian, it 
should be clear, due to the obvious presence of oscillatory terms, that the system does 
not satisfy the adiabatic approximation even if Eq.([T]) holds. Contrarily, when no such 
analytical expression is known, the inconsistency may be more relevant. This case occurs 
when the Hamiltonian is known only through experimental information on its energy levels 
and eigenstates. 

In conclusion, the present work eliminates the mystery that has surrounded the adiabatic 
theorem since the publication of Ref . [5j . Although it is expected that only a few systems 
will be affected by the inconsistency the condition Eqs. fl22l) . flMl) must supplement the 
quantitative conditions Eq. ([T]) whenever the validity of the adiabatic theorem and/or 
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adiabatic approximation are investigated. 
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